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Abstract — We demonstrate a technique for shaping current 
inputs for the direct modulation of a semiconductor laser for 
digital communication. The introduction of shaped current 
inputs allows for the suppression of relaxation oscillations and 
the avoidance of dynamical memory in the physical laser device, 
i.e., the output will not be influenced by previously communicated 
information. On the example of time-optimized bits, the possible 
performance enhancement for high data rate communications is 
shown numerically. 

Index Terms — Semiconductor lasers, Intersymbol interfer- 
ence, Time optimal control, Optimal control, Communication 
system nonlinearities, Nonlinear optics, Nonlinear systems 



I. Introduction 

SEMICONDUCTOR lasers typically convey information 
for optical communication using either "external" or "di- 
rect" modulation. In the first scheme, the laser is maintained 
in a constant light-emitting state and some sort of exter- 
nal "shutter," typically a Mach-Zehnder type electro-optic 
modulator, modulates the output intensity according to an 
externally applied voltage. The advantage of this technique 
is that it allows high bit rates (> 10 Gbit/s) and there 
is no interaction between the information being applied to 
the modulator and any internal dynamics of the laser. The 
disadvantage is that costly external apparatus are necessary. 
Direct modulation involves changing the current input of the 
laser, i.e., modulating it around the bias current, in order to 
produce time-dependent output in optical intensity. 

In terms of modeling, the use of direct modulation means 
that the laser becomes a driven nonlinear system with po- 
tentially complicated dynamics. Indeed, it is well known that 
even simple high-frequency periodic modulation of sufficiently 
large amplitude can lead to nonlinear phenomena such as 
period doubling cascades, period tripling, and chaos [l]-[8]. 
It is therefore not surprising that attempting communication 
with current pulses that occur rapidly in succession, quickly 
changing among discrete current levels to transmit randomly 
interleaved zero and one bits, will result in a complex output 
time series of intensity. As a result the bits may cease to be 
decodable at the receiver and there is a considerable memory 
effect [9]. This effect is called "inter-symbol interference" 
in the communications literature because the sequence of 
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information symbols at previous times may harm decoding 
at the present time. Additionally, communication using high- 
speed direct modulation is limited by transient oscillations, 
known as relaxation oscillations, that result from the interplay 
between the optical field and the carrier density [10], [11]. 
Such oscillations cause distortions of the output light pulse 
and broaden the signal's optical spectrum (increased chirp), 
thus contributing to an increased bit-error rate. 

In an effort to overcome the problems associated with 
direct modulation many techniques for chirp reduction and 
suppression of relaxation oscillations have been proposed, 
including the modification of the physical laser device [11]- 
[13], the use of external optical feedback [14] and external 
electronic circuits [15], as well as the shaping of the current 
pulses that are used to encode logical bits into the optical 
output of the laser [16]-[20]. 

We approach the issue of inter-symbol interference and 
relaxation oscillations from the point of view of nonlinear 
dynamics. A basic single-mode semiconductor laser modeled 
by the rate equations is a two degree of freedom system, 
with ordinary differential equations for photon density and 
carrier density. The cause of inter-symbol interference is the 
fact that even if the observable intensity returns to its nominal 
stationary condition after a current wave form is applied, the 
unobservable carrier density may not have, and thus there 
is dynamical coupling from past to future. Our feed-forward 
solution is to use appropriately shaped driving currents to 
control the state-space trajectory of the laser, and we show 
how to obtain these waveforms given a dynamical model for 
the laser. 

The main goal of this paper is to demonstrate a tech- 
nique that can eliminate relaxation oscillations and inter- 
symbol interference, without having to physically change the 
laser diode. This allows enhanced data rates and improved 
communication performance as is shown numerically for a 
communication scheme using time-optimized non-return-to- 
zero (NRZ) bits. This improvement comes at the cost of 
more complex input electronics, which would be necessary 
to shape the input current in experiments. A distinguishing 
feature of our approach is the use of shaped, continuous, drive 
currents with finite rise times. This contrasts with other recent 
research [20], where relaxation oscillations are suppressed 
using piece-wise constant currents, which would possibly 
present difficult discontinuities. 

We will first present the general idea of the shaping method 
in Sec. |n] and then show particular solutions for the specific 
case of optimizing the transition time for NRZ bits in Sec.lIIII 
Using these solutions, we estimate the bit error rate (BER) 
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performance for the case of a simple Gaussian white noise 
channel to show the possible performance enhancement. 

II. Current Shaping 

The semiconductor laser is presumed to be single mode. For 
this case the laser dynamics can be described by a driven two 
dimensional system of ordinary differential equations for the 
carrier density N(t) and the photon density S(t), the latter 
being proportional to output intensity: 

= - lc S + TG(N,S) S, 

(1) 
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dt ed 
The constant Jo is the bias current, J(t) the modulation around 
the bias, G(N(t), S(t)) the optical gain coefficient including 
nonlinear effects, T the confinement factor, d the active layer 
thickness of the laser, j c the photon decay rate, and 7„ the 
spontaneous carrier decay rate. 

For the purpose of numerical calculations and analytical in- 
vestigations it is useful to work with dimensionless quantities. 
It is convenient to non-dimensionalize using the fixed point 
(So, -Wo) °f the laser in the absence of modulation ( J(t) = 0), 
that is, of the laser under stationary continuous-wave operating 
conditions (dS/dt = dN/dt = 0). We expand the nonlinear 
G(N(t), S(t)) around this point to first order 



G(N, S) = G + G n (N - N ) + G P (S - S ). 



(2) 



Liu and Simpson [21] show how to experimentally estimate 
G n = dG/dN > and G p = dG/dS < 0. Note that the fixed 
point conditions yield the relations TGq = j c and Jq /ed — 
j s N = G S . 

These dimensionless dynamical equations are independent 
of how far above threshold the laser is operated except via 
indirect influence on the empirical differential gain parameters 
G n and G p , which depend on the expansion point. The dimen- 
sionless photon density s(t) is defined via s(t) = S(t)/So, the 
carrier density n(t) via n(t) = N(i)/No — 1, and the current 
via J + Jm = {Jo + J(t))/{j s N ed) — l. [ The dimensionless 
gain g = G/Go is 
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where 7„ = G„So and j p = —TG p Sq. Using the natural 
relaxation oscillation angular-frequency ujr — y'jdn 
the dimensionless quantity corresponding to time r 
be formed. With this the model equations are: 

ds _ 7c 
dr ur 
dn _ 7 S 
dr u R 

Table U lists the dimensionless quantities and the parameters 
we use, the values of which are based on measurements of J.M. 
Liu's group [22]. For the purpose of digital communication 

'Roughly speaking, the dimensionless total current 1 + J + i s m 
units of the threshold current. This is strictly true for simple gain functions 
that imply clamping of the carrier density to it's threshold value for currents 
above threshold. 



TABLE I 

Dynamical Variables, definitions, and numerical values of 
parameters. 



Symbol 


Value 


Description 


s{t) 




dimensionless photon density 


n(r) 




dimensionless carrier density 


g(n,s) 




dimensionless gain coefficient 


Jm(t) 




dimensionless modulation current 


J 


2/3 


bias current at fixed point 




1.458 xlO 9 ^- 1 


spontaneous carrier decay rate 


7n/J 


2.0 xl0 9 s" 1 


gain variation with carrier density 


lv/J 


3.6 xlO 9 ^" 1 


gain variation with photon density 


7c 


3.6 xlO 11 ^- 1 


photon decay rate 




V7c7n + 7s 7p ~ 3^: 


relaxation oscillation frequency 



we are interested in the relationship between the information 
bearing modulation current and the output of the laser. The 
key insight is to note that, given a functional form for the time 
series s(t), the corresponding input current can be computed 
just by differentiation. In contrast, given the functional form of 
the current, it is necessary to integrate the nonlinear equations 
of motion (|4j to obtain the potentially chaotic resulting laser 
output. 

Given a particular model for the laser, there is often some 
differential operator taking a given s(r) to the other variables 
O s \s{t)\ = Jm{ t ) an d another one O n \s(r)\ = n(r) such 
that the state vector [JM{ T ), n { T ), S ( T )] lS a solution to 
the equations of motion. This solution though may not be 
physically allowable (such as having negative total current) 
or dynamically stable. 

To derive the driving current as function of the output 
waveform we introduce, for computational ease, s(t) = e ytyT \ 
solve for n(r) using the ds/dr equation, differentiate, and 
rearrange the dn/dr equation to solve for Jm{ t ), yielding 



Jm 



Icln 



Ip + ln 



— )y+(e y 



i: 



(5) 



This equation relates a given waveform s(r), via y, to the 
modulation current that caused it. Any desired output wave- 
form that corresponds to a physically allowed input can thus 
be produced by driving the laser with the according shaped 
current Jm{ t )- Not all functions s(t) will however result in 
realistic Jm{t)- One physical requirement is a non-negative 
total driving current at all times, 1 + J + Jm > 0, as 
it is not possible to extract electrical energy back from the 
carrier density. Practical considerations generally limit the 
peak current that can be applied as well. 

Note that l|5} includes the full nonlinear dynamics of the 
laser as modeled by the rate equations Q. In contrast to any 
theory based on the linearization of Q, results obtained using 
(|3 apply to large driving amplitudes. The specific linearized 
parameterization of gain we used restricts the validity 
of because it does not take into account gain saturation 
effects that become important for very large intracavity photon 
densities s. However, it is equally straightforward to derive 
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Time (in units of l/co R ) 



Fig. 1 . We show how a square-like total driving current 1 + J + Jj^ (a) 
gives rise to relaxation oscillations of the laser photon density s(t) (b). The 
bit time is Tj,; t = 8 which corresponds to 360 ps for the model-parameters 
shown in Tab. Qlu;^ 1 ~ 45 ps). The lower and upper value for the driving 
current corresponds to fixed points with s = 0.2 and s = 2.0, respectively. 



formulas analogous to 10 when gain functions which model 
saturation are used. We chose to work with gain (0 because 
it is simple. For the results presented in the rest of the paper 
gain saturation is not important, and, the parameters of 
may be estimated experimentally. 

III. Application to Communication 

We consider a "direct" modulation setup: the laser's driving 
current J(t) is modulated to encode information. The laser 
light travels through the communication channel, and at the 
receiver, the intensity of the arriving signal is measured. One 
then attempts to extract the information from this signal, in 
the simplest case by a single thresholding to decode binary 
information. 

Figure n shows the laser dynamics driven with an unshaped 
current waveform (NRZ bits). After a transition from low 
current to high current, there are clear relaxation oscillations 
of the laser output, that is, the output "rings" at the char- 
acteristic relaxation oscillation frequency. In the particular 
situation here, the bit time would be sufficiently large that a 
conventional system could decode the information successfully 
by integrating energy over a bit interval. A further increase 
in bit rate, however, will result in significant inter-symbol 
interference (see Fig. |3jl-f) 

The state-space trajectory depends on its initial conditions 
s° and n° as well as on the driving signal Jm{ t )- If, at the 
end of each bit time Tuu the system returns to a particular 
point in state-space, e.g. s(T^ t ) = s° and n(Xfo t ) = n°, 
then repeating the same driving signal Jm{ t ) will result 
in a repetition of the state-space trajectory. If, on the other 
hand, the system trajectory did not reach (s°,n°) after the 
first bit, then repeating the driving signal will result in a 
different laser output, whose particular details depend on the 
previous drive signal. This dynamic coupling causes inter- 
symbol interference. 

The main physical issue is to return the unobserved carrier 
density to a specific state as well as the photon density. We 



use to design driving signals which return the system 
trajectory to a known state at the end of each bit. Since 
the adverse effects of dynamic coupling on communication 
performance are especially relevant for high bit rates, we will 
try to maximize the bit rate (minimizing the bit time) for a 
given physical laser by optimizing driving signals which cause 
no dynamical inter-symbol interference. 

A. Finding the time-optimal NRZ bits 

NRZ bits encode 1-bits (0-bits) as high (low) light power, 
but in contrast to return-to-zero pulsing, the intensity stays 
high when two or more consecutive 1 -bits are transmitted, and 
similarly for consecutive zero bits. The dynamical picture is a 
system that can be driven to switch between two fixed points 
in state-space, one corresponding to low output power, the 
"off" state, and the other corresponding to high output power, 
the "on" state. The driving waveform for a bit identical to its 
predecessor is trivial — a constant to stay at the particular fixed 
point. We solve for driving currents that fix the system along 
a state-space trajectory, smoothly connecting the "off" fixed 
point to the "on" fixed point (up transition) and another one 
for the "on" to "off" transition (down transition). We require 
both the photon density and the carrier density come to a 
known fixed point at the end of each bit, thus there will be no 
dynamic coupling from one bit to the next and therefore no 
inter-symbol interference. 

We use (0, a parameterized set of photon density wave 
forms, and an optimization algorithm to search for waveforms 
which optimize some other desirable communication criterion 
under the condition that they correctly connect the two states. 
In particular, we have searched for the minimum transition 
time T, and thus highest bit rate, given a desired intensity 
separation between high and low. 

To be specific, let us fix the "off" intensity to 20% of the 
fixed point intensity around which we normalized, s s = 0.2, 
and the "on" state to twice the fixed point value, s on = 2.0. 
For the up (down) transition we need to search the space of 
all smooth waveforms s(r) with s(0) = s ff and s(T) = s on 
(s(0) = s on and s(T) = s s) such that the transition time T 
is minimized and such that the corresponding current Jm (t) 
meets certain constraints. In this example, we require that < 
1 + J + Jm < 6, which means that the total input current 
into the laser is positive and that the largest pump current we 
admit is roughly six times the threshold current, a physically 
reasonable requirement. We furthermore require that Jm{t) 
has zero first derivatives at r = and r = T. 

The success of the optimization over all possible waveforms 
depends on the choice of the right set of basis functions, and 
satisfying all the boundary conditions. We found that, although 
an ordinary Fourier base does work to a certain extent, it 
is numerically advantageous to pick base functions which 
automatically incorporate many of the constraints, avoiding 
needing as many explicit constraint equations in the opti- 
mization procedure. This notion led us to use the following 
functions to parameterize ,s = e v : 

-, 4 

^ a k sin ( 2nf 
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Fig. 2. The figure shows driving current waveforms, i.e. 1 + J + Jm. M> 
and in the insets the corresponding optical laser output in terms of the photon 
density s(t). Panels a), b), and c) show the transition from the "off" to the 
"on" state. Using two sin-terms [N=2 in |6|] in the optimization procedure 
for the driving current yields a), whereas the use of 20 base functions (N=20) 
results in b). For comparison, we display in c) the photon density waveform 
and discontinuous current computed using optimal control theory. Panels d), 
e), and f) are results for the "on" to "off" transition, where two base functions 
are used in d), 20 base functions in e), and f) is the optimal control theory 
result. 

where p(-) is a (7th degree) polynomial with p(0) = In [s(0)], 
p(l) = In [s(T)], and with zero first, second, and third deriva- 
tive at t = and r = T, which ensures zero first derivatives 
of Jm at these times. f(x) = (a + l)x/(aox + 1) is a 
warp function introducing asymmetry with the effect that short 
transition times can be obtained using fewer base functions. 
The only constraint that has to be enforced in the optimization 
is the requirement < 1 + J + Jm < 6 (see Table |H}. 
This choice of base functions was rather arbitrary, justified 
principally by numerical ease in automatically incorporating 
boundary conditions and empirically giving reasonably suc- 
cessful results, not by any strong theoretical considerations. 
We leave the question of how to choose mathematically 
optimal, or experimentally desirable, base functions for future 
investigations. Our main goal is to demonstrate the general 
concept; there is no barrier to incorporating other forms of 
base functions and other models for different lasers in this 
method. 

Figure I2 shows numerical results from using base functions 
(|6j and the NPSOL software package [23], [24], which allows 
nonlinear constrained optimization. Comparing the currents 
for the "off" to "on" transition (Fig. [2^, b, and c) to the 
corresponding photon densities s(i) (insets) shows that the 
shaped current waveform that is obtained as result of the 
optimization is quite dissimilar from the output intensity 
waveform. 

In particular, an increase of the laser's output intensity is not 
achieved by just ramping up the current from the low value 
corresponding to the "off" steady state to the higher value 
corresponding to the "on" state. Instead the pump current is 
reduced after an initial surge to allow the carrier density to 



TABLE II 

Summary of parameters, constraints, and results for the 
optimized nrz bits 



Fixed parameters in the optimization 


0- bit fixed point 

1- bit fixed point 

minimal allowed pump power 
maximal allowed pump power 


sout = 0.2 

*lbit = 2 -0 
1 + J + JM min =0.0 
1 + J + J M max = 6.0 


Optimization constraints and how they are met 


fixed s(0) (either s on or s„jf) 
fixed s(T) (either s D ff or s„n) 

f t JM(T) = 

<1 + J + J m (t) < 6 


2/(0) = p(0) = ln(s(0)) 
y(T) = p(l) = ln(«(T)) 
y(0) = y( ) = y(0) = 
y(T) = y(T) = y(T) = 
active constraint 


Results of optimization 




N=2 


N=20 


Theory 


up transition 
down transition 


T = 1.80 
T = 1.98 


T = 1.57 
T = 1.71 


Topt = 1.42 
T op t = 1.53 



relax to the steady state value. 

In Fig. |2ji we show the current and photon density that 
were obtained with just two sinusoidal terms (N=2 in (J6J). 
The minimal transition time that can be obtained decreases 
when the number of base functions is increased, e.g. tenfold 
(N=20) as shown in Fig.|2j). The inclusion of higher frequency 
sine functions results in a current with a much faster rise time. 
Just by the choice of the number of base functions one can 
therefore approximately take into account limitations on the 
possible rise-times of electronics involved. Thus, we did not 
explicitly add this limitation as an additional constraint in the 
optimization procedure, but such a constraint is reasonable to 
use. 

The waveforms found also converge to the theoretical best 
limit (Fig. |2fc), obtained by applying results from optimal 
control theory. If one drops the requirement of finite rise times 
for the drive current, that is, if discontinuous functions are 
admitted, then the so called Minimum Principle [25] can be 
applied. It tells us that a time-optimal state-space trajectory 
is achieved by driving the laser with a constant current at 
either the maximal allowed or minimal possible value (see as 
well [20]). The point in state-space where to switch from 
one value to the other can be found by integration of the 
Hamiltonian corresponding to @ (for details see [26]). The 
resulting minimum transition time is a lower bound benchmark 
for the transition times that can be achieved with realistic input 
currents. It also is the lower bound for the bit period T^u = T 
for which communication without laser induced inter-symbol 
interference is possible. Figure [2] shows that the results we 
obtain with our technique converge both in the shape of the 
waveform and in the value of the transition time to the optimal- 
control theory result as more base functions are included in 
the optimization routine. 

The down transition driving the laser from the "on" power 
level to the "off" power level (Fig. |2jl, e, and f) we 
design exactly as before except that values used for the 
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boundary conditions s(0) and s(T) are reversed from the low 
to high transition. The minimal transition times are displayed 
in table |H] where we show the results using two and twenty 
base functions as well as using optimal control theory. It 
is important to note that the down transition is not just a 
time reversed version of the up transition, and that both the 
transition times and the waveforms are different and have 
to be designed separately. Most communication schemes are 
evenly clocked, the time necessary to transmit a symbol is 
independent of the particular information. Yet our results show 
that the fundamental dynamics of the laser may mean that 
approaching the optimal limits might imply using unequal 
symbol times. We offer this as an idea to explore but we do 
not pursue it further. 

B. Performance 

Using shaped input currents will result in a better BER 
performance for high bit rates because we avoid errors caused 
by inter-symbol interference due to the internal dynamics 
of the laser. To show the possible improved communication 
performance we compare NRZ bits obtained using the up and 
down transitions that were designed by the presented technique 
using two base functions (N=2 in (|6j) to unshaped current 
modulation. 

A random sequence of bits is to be transmitted in this 
numerical experiment. The bit sequence determines the shaped 
current, shown in Fig|3ji, which enters as driving term the laser 
rate equations (0}. The output s(r) is obtained by integration 
of @ and Fig|3j5 shows that s(r) switches, as was the design 
goal, between the "on" and "off" state without any relaxation 
oscillations. Transmitting 50 bits results in the state-space 
trajectory presented in Fig|3j; (thick line). The two fixed points 
(square and diamond) and the optimal control trajectory (dash- 
dotted line) obtained by using non-continuous currents are also 
shown. 

The bit time for all the panels in Fig[3]was set to Tut = 2.0. 
We extend each transition to this time by appending the ap- 
propriate interval of fixed-point current. Since we normalized 
time using lur = 2nf r this implies that we communicate using 
NRZ bits that last about one third of the relaxation oscillation 
period l/f r of the fixed point (So, No). 

We compare the bit errors using shaped currents to square- 
like current modulation where, basically, the current follows 
the NRZ bits but with finite rise times, to be slightly more 
realistic. The up and down transitions of the square-like 
driving current shown in Fig|3jl was achieved using 'tanh' 
functions. The maximum slope of the current equals the 
maximum slope of the shaped current shown in Fig|3^. Since 
the bit time is short compared to the period of the laser's 
intrinsic relaxation oscillations, the state-space trajectory due 
to the modulation of 50 NRZ-bits is complicated and exhibits 
large excursions (see Fig|3f). Since the system does not reach 
either one of the two fixed points at the end of a bit-period, 
there is considerable dynamic coupling from past to future 
and hence inter-symbol interference. In Fig|3^ we show the 
output corresponding to the drive current above (Fig|5Jl). It is 
apparent that even without channel noise the receiver will not 
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Fig. 3. Comparison between shaped current modulation (a-c) and square- 
like currents (d-f). Shaped drive currents (a) result in an optical laser output (b) 
free of relaxation oscillations. For visual guidance, bit times are indicated by 
dotted lines. Note that the currents were designed to optimize the transition 
between the two logical states ("on" and "off" ) and consequently, in the 
output, the assigned bits appear shifted by half a bit time. The bit assignment 
is shown between panels a) and b) together with the corresponding bit energies 
(see text). Square-like currents (d) result in an output (e) corrupted by strong 
inter-symbol interference. Panels c) and f) are the corresponding state-space 
plots (see text for further details). 



be able to infer without error the bit sequence that was to be 
transmitted. 

To quantify the possible improved communication perfor- 
mance we calculated bit error rates (BER). For illustration, a 
simple additive white Gaussian noise corrupted a dispersion- 
less channel. The overall conclusions drawn will apply to more 
realistic channel models, since the performance enhancement 
is due to avoiding inter-symbol interference in the laser and 
should be largely independent of the characteristics of the 
communication channel. 

Figure [4] shows in the main panel the SNR versus the BER 
for high bit rate (Tut = 2) square-like current modulation 
as squares and high bit rate shaped current modulation as 
circles. Low bit rate (Tut = 8) square-like current modulation 
(squares) and low bit rate shaped current modulation (circles) 
is presented in the inset of Fig. |4] These are results of nu- 
merical experiments where one million random bits generated 
a corresponding input current (either shaped or square-like) 
to @, which was integrated. To the resulting laser output an 
appropriate amount of channel noise was added. The bits were 
decoded at the receiver by integrating the received power over 
each bit time (assuming a perfect clock) and by distinguishing 
zero and one bits through simple thresholding. This way an 
estimate of the BER was obtained for each noise level. 

For the case of square-like current modulation the threshold 
used for bit detection was determined numerically. It was 
optimized for each noise level separately. Furthermore, due to 
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SNR in dB 

Fig. 4. The signal-to-noise ratio (SNR) versus the bit-error-rate (BER) 
is presented for the laser being driven by both shaped currents (circles and 
lines) and square-like currents (squares). The main panel shows high bit rate 
(T&it = 2) and the inset lower bit rate (T^ it = 8) communication. The 
symbols represent results obtained by numerical integration and the lines are 
theoretical curves. We display as triangles results where stimulated emission 
noise is simulated. 



was used. 

Calculating the theoretically expected BER using the opti- 
mal threshold value E?£ (a), which can be derived from @, 
and which is a function of the noise strength, results in a curve 
that to the eye is indistinguishable from the one displayed in 
FigE For this reason we used the simpler E t n given by (jSJl 
in the numerical experiment determining the BER when using 
shaped modulation currents. 

The following conclusion can be drawn from the results pre- 
sented in Fig.|4] As expected, using shaped current modulation 
there is a very good agreement between the theoretical and 
numerical result for both high and low bit rates, showing that 
all errors are due to Gaussian channel noise. In contrast, Fig.0] 
shows that high bit rates (Tut = 2) for square-like current 
modulation lead essentially to a communication breakdown. 
Even for very small noise levels (high SNR) on average about 
14 bits per 100 transmitted bits are wrongly decoded by the 
receiver. These errors are caused by inter-symbol interference 
due to the laser dynamics. For lower bit rates, that is bit rates 
smaller than w r /2n, the bit errors attributable to the laser 
dynamics become more and more negligible, as is indicated 
by the inset of Fig. 0] which shows that for lower bit rates 
the performance enhancement due to the use of shaped input 
currents is minor. 



the complex dynamics of s(r), it is not immediately obvious 
when to start the receiver clock that determines the beginning 
and end of a bit (see Fig. [3^). Therefore the optimal starting 
time was numerically estimated as well. 

For a shaped current input, a simple theoretical expression 
for the expected BER can be derived. Consider transmission 
without channel noise first. Without any noise the receiver 
detects four different bit energies E l (i £ {1,2,3,4}) sig- 
nifying one bits. 2 The largest energy E 1 corresponds to a 
one bit being both preceded and succeeded by a one bit, 
the smallest energy E A corresponds to a one bit being both 
preceded and succeeded by a zero bit, whereas E 2 and E 3 
correspond to one bits that contain an up or down transition 
(see Fig. Similarly, there are four energy levels (E^ 
with j G {5, 6, 7, 8}) representing zero bits. The expected bit 
error rate is obtained by evaluating the probability that the 
energy that results after adding white Gaussian noise £ with 
variance a 1 to a zero bit energy level exceeds the threshold 
E t h and by considering the probability that E l + £ < E t h for 
i € {1,2, 3, 4}. Doing this the theoretically expected BER can 
be written in terms of the complementary error function: 

(7) 
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The resulting curves are plotted as lines in Fig 
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2 For the purpose of the numerical experiment we use E = ^ htt s(t)dt 
for the bit energies because the photon density s(t) is proportional to the 
laser output power. 



C. Noise and Parameter Mismatch 

Semiconductor laser diodes are inherently noisy. To study 
the effects of spontaneous emission noise on the proposed 
communication scheme we conduct numerical experiments, 
where rate equations that include appropriate stochastic terms 
(Gaussian white noise) are driven by the shaped current, 
which, as discussed above, was calculated based on the deter- 
ministic model (|4}. The rate equation including spontaneous 
emission terms and a description of the integrator can be found 
in [22], [26]. 

For realistic values of the noise variance (same as in [22]) 
the performance advantage achieved by modulating the laser 
with shaped input currents persists. Although noise causes the 
state-space trajectory to deviate from the ideal trajectory, the 
deviations remained for the bit sequences and noise realization 
of our numerical experiments small enough to cause only 
slight performance degradations. The BER performance calcu- 
lated in one such numerical experiment is shown as triangles 
in the main panel of Fig. 0] It is close to the theoretical curve 
calculated for the noise free system and far better than square- 
like modulation of a laser model free of spontaneous emission 
noise. 

The shaped current waveform that is used to achieve high 
bit rate communications without inter-symbol interference was 
calculated based on the model @, the parameters of which 
will be known only approximately for a given laser diode. 
It is therefore important to consider the effects of parameter 
mismatch between the model used to generate the shaped 
driving current and the model used to numerically evaluate 
the BER performance. 

Figure |5] shows the degradation of the BER performance for 
differential gain mismatches, defined as A = (7,/ — j n ) /j n '- 
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SNR in dB 

Fig. 5. BER performances for mismatches A = ? n , where -y n ' is the 
estimated and 7 ?l the actual differential gain. The inset shows a time trace of 
the photon density s for A = —10% at the output of the laser. The bits are 
clearly distinguishable. 

Here j n ' denotes the estimated differential gain, the one used 
in the generation of the shaped drive current, and 7„ is the 
actual parameter, the one characterizing the system driven by 
the shaped currents. Overestimating of the differential gain 
(A > 0) severely degrades performance for mismatches | A| > 
5%, whereas underestimation (A < 0) has a far less noticeable 
effect. 

The reason is that overestimating the differential gain im- 
plies an overestimation of lor, i.e. the response time of the 
actual laser is slower than was assumed. The laser will have 
slower rise and fall times and consequently a lower bit energy 
for a one bit that is both preceded and succeeded by a zero bit. 
Although the mean bit energy does not change significantly, 
the energy difference between the maximum energy for a zero 
bit (max j(E [5 , 8 ](£ J )) and the minimum energy for a one bit 
(min^R M(E % y) is decreased resulting in a higher BER. 

A similar asymmetric performance degradation holds for 
mismatches in the photon decay parameter j c . Performance is 
most sensitive to mismatches of j c and 7„ and least sensitive 
to mismatches in 7 p . 

The unexpected asymmetry of the performance degradation 
with estimation may suggests that, in practice, it may be 
wiser to solve for state-space trajectories using slightly smaller 
values of 7 Ci „ than the best available experimental estimations 
of the parameters. This may reduce the risk to performance 
from fluctuations in 7 C „. 

IV. Discussion 

We demonstrate a technique to control the laser's nonlinear 
internal degrees of freedom to move along a specific state- 
space trajectory by designing current input waveforms. This 
can be used to optimize the optical intensity waveforms 
emitted to fit certain desired characteristics. 




Time (in units of 1/co ) 



Fig. 6. Evolution of the photon density s(t) in time for initial conditions 
chosen randomly in the neighborhood of the correct initial state (a subset 
of ten realizations is displayed, no noise and no mismatch). Inset: Initial 
conditions in state-space (circles) and where the trajectories in state-space 
end up (squares) after driving currents for 49 random bits were applied. 

On the example of time-optimized NRZ bits, the possible 
BER performance enhancement due to the elimination of inter- 
symbol interference is demonstrated. We furthermore show 
that spontaneous emission noise and parameter mismatches 
can be tolerated suggesting that this scheme can be used in 
practice to enhance the communication data rate when directly 
modulated lasers are used as transmitters. 

Our numerical experiments show that including spontaneous 
emission noise and small but finite parameter mismatches 
leads to trajectories in state-space that are close to the ideal 
one. This suggest that the trajectory we obtain in the opti- 
mization is robust with respect to perturbations. In numerical 
investigations it appears that the state-space trajectory is also 
stable in the sense that, starting from slightly perturbed initial 
conditions, trajectories converge to the ideal one, when the 
system is driven by the same bit sequence. The result of 
one such numerical experiment is shown in Fig. |6] where 
initial relaxation oscillations soon disappear as the trajectories 
converge. 

Theoretically, an open dynamical-systems question is 
whether there is some kind of easily computable generalized 
stability property that would quantify whether the controlled 
system is stable even under arbitrary message driving. It 
would be desirable to build this property into the optimization 
target to ensure that the solutions found are also stable. In 
our model it was fortuitous that all solutions were stable 
with the parameters we used, but this was not guaranteed 
by the method, and at present the best check is empirical 
integration under perturbations. We also speculate that the 
relative advantage of dynamically-aware current-shaping over 
conventional direct modulation could be significantly larger 
if, instead of transmitting binary symbols (two energy levels), 
the system was controlled to a higher alphabet. The net 
result of our pulse shaping is to radically reduce that part of 
variance which is due to dynamics in the optical intensities 
of the different received symbols. If the fluctuations from 
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inter-symbol dynamics are eliminated, leaving only detector 
noise and spontaneous emission, then it is conceivable that 
one could transmit with more than two distinct optical levels 
for the symbols, and still reliably decode them. This would 
improve useful data rate versus binary transmission even 
without increasing the rate of symbols being transmitted. 

The key technological issue is whether the proposed direct- 
modulation scheme using shaped currents can be implemented 
at low cost. Existing technology has been developed for clock 
equalization, generating fast fronts at gigahertz speeds which 
could subsequently be shaped using delay line techniques [27]. 
Thus, current shaping is in principle electrically achievable. 
We speculate that the laser and the current shaping electronics 
could be integrated on one chip so that the overall cost should 
be comparable to that of traditional single-mode-lasers. We 
hope that the potential of current-shaping to mitigate problems 
that have plagued the direct-modulation approach to high- 
speed communication will stimulate further research and lead 
to practical applications. 

References 

[1] C. H. Lee, T. H. Yoon, and S. Y. Shin, "Period doubling and chaos 
in a directly modulated laser diode," Appl.Phys.Lett., vol. 46, no.l, pp. 
95-97, Jan. 1985 

[2] Y. C. Chen, H. G. Winful, and J. M. Liu, "Subharmonic bifurcations 
and irregular pulsing behaviour of modulated semiconductor lasers," 
Appl.Phys.Lett., vol. 47, no.3, pp. 208-210, Aug. 1985 

[3] M. Tang and S. Wang, "Simulation studies of bifurcation and chaos 
in semiconductor lasers," Appl.Phys.Lett., vol. 48, no. 14, pp. 900-902, 
April 1986 

[4] G. P. Agrawal, "Effect of gain nonlinearities on period doubling and 

chaos in directly modulated semiconductor lasers", Appl.Phys.Lett., vol. 

49, no.16, pp. 1013-1015, Oct. 1986 
[5] W. Lauterborn and R. Steinhoff, "Bifurcation structure of a laser with 

pump modulation," /. Opt. Soc. Am. B, vol. 5, no. 5, pp. 1097-1104, 

May 1988 

[6] H. Lamela, G. Carpintero, and F. J. Mancebo, "Period tripling and chaos 
in the dynamic behavior of directly modulated diode lasers," IEEE J. 
Quantum Electron., vol. QE-34, no. 10, pp. 1797-1801, Oct. 1998 

[7] K. A. Corbett and M. W. Hamilton, "Comparison of the bifurcation 
scenarios predicted by the single-mode and multimode semiconductor 
rate equations," Phys. Rev. E, vol. 62, no. 5, pp. 6487-6495, Nov. 2000 

[8] H. F. Liu and W. F. Ngai, "Nonlinear dynamics of a directly modulated 
1.55 p,m InGaAsP distributed feedback semiconductor laser," IEEE J. 
Quantum Electron., vol. QE-29, no. 6, pp. 1668-1675, June 1993 

[9] P. Colet, C. R. Mirasso, and M. S. Miguel, "Memory Diagram of Single- 
Mode Semiconductor Lasers", IEEE J. Quantum Electron., vol. QE-29, 
pp. 1624-1630, June 1993 
[10] K. Petermann, Laser Diode Modulation and Noise (Kluwer, Tokyo, 
1988) 

[11] K. Y. Lau and A. Yariv, "Ultra-High Speed Semiconductor Lasers", 

IEEE J. Quantum Electron., vol. QE-21, pp. 121-138, Feb. 1985 
[12] R. S. Tucker, "High-Speed Modulation of Semiconductor Lasers", /. 

Lightwave Techno!., vol. LT-3, pp. 1180-1192, Dec. 1985 
[13] S. Riyopoulos, "Elimination of transient vertical-cavity surface-emitting 

laser oscillations using photoactive feedback", Appl. Phys. Lett, vol. 75, 

pp. 3057-3059, Nov. 1999 
[14] C. R. Mirasso, E. Hernandez-Garcia, J. Dellunde, M. C. Torrent, and J. 

M. Sancho, "Analytical calculations of switch-on time and timing jitter 

in diode lasers subjected to optical feedback and external light injection", 

Opt. Comm., vol. 115, pp. 523-527, April 1995 
[15] Y. Suematsu and T. H. Hong "Suppression of Relaxation Oscillation 

in Light Output of Injection by Electrical Resonance Circuit", IEEE J. 

Quantum Electron., vol. QE-13, pp. 756-762, Sept. 1977 
[16] M. Danielsen, "A Theoretical Analysis for Gigabit/Second Pulse Code 

Modulation of Semiconductor Lasers", IEEE J. Quantum Electron. , vol. 

QE-12, pp. 657-660, Nov. 1976 



[17] P. Torphammar, R. Tell, H. Eklund, and A. R. Johnston, " Minimizing 
Pattern Effects in Semiconductor Lasers at High Rate Pulse Modula- 
tion", IEEE J. Quantum Electron., vol. QE-15, pp. 1271-1276, Nov. 
1979 

[18] R. Olshansky and D. Fye, "Reduction of dynamics linewidth in single- 
frequency semiconductor lasers", Electron. Lett., vol. 20, pp. 928-929, 
Oct. 1984 

[19] L. Bickers and L. D. Westbrook, "Reduction of transient laser chirp in 

1.5 ^m DFB lasers by shaping the modulation pulse", IEE Proceedings, 

Pt. J, vol. 133, pp. 155-162, April 1986 
[20] N. Dokhane and G. L. Lippi, "Improved direct modulation technique 

for faster switching of diode laser", IEE Proc.-Optoelectron. , vol. 149, 

pp.7- 16, Feb. 2002 

[21] J. M. Liu and T. B. Simpson, "Four-Wave Mixing and Optical Modula- 
tion in a Semiconductor Laser", IEEE J. Quantum Electron., vol. QE-30, 
pp. 957-965, April 1994 

[22] H. D. I. Abarbanel, M. Kennel, L. Illing, S. Tang, H.F. Chen and J.M. 
Liu, "Synchronization and Communication Using Semiconductor Lasers 
With Optoelectronic Feedback", IEEE J. Quantum Electron. , vol. QE- 
37, pp. 1301-1311, Oct. 2001 

[23] P.E. Gill, W. Murray, M.A. Saunders, and M.H. Wright, "Some theoreti- 
cal properties of an augmented Lagrangian merit function", Advances in 
Optimization and Parallel Computing, P.M. Pardalos, ed., North Holland, 
1992, pp. 101-128. 

[24] P.E. Gill, W. Murray, M.A. Saunders, and M.H. Wright, "User's Guide 
for NPSOL (Version 4.0): A Fortran package for nonlinear program- 
ming", Report SOL 86-2, Department of Operations Research, Stanford 
University, Stanford, CA, 1986 

[25] M. Athans and P. L. Flab, Optimal Control: An Introduction to the 
Theory and Its Applications, McGraw-Hill Book Company, 1966 

[26] L. Illing, Chaos synchronization and communication in semiconductor 
lasers. Thesis (Ph.D.), University of California, San Diego, 2002 

[27] W. J. Dally, J. Poulton, "Transmitter equalization for 4-Gbps signaling", 
IEEE Micro, vol. 17, pp. 48-56, Jan. 1997 




Lucas Illing was born in Halle, Germany, in 1973. 
He received his Vordiplom in 1996 from the Hum- 
boldt University in Berlin, Germany and in 2002 the 
Ph.D. in Physics from the University of California, 
San Diego. He is a Research Assistant at Duke 
University. His research is primarily concerned with 
chaos synchronization, control, and communication 
using chaos. 



Matthew B. Kennel was born in Santa Monica, California in 1968. He 
received the A.B. degree in physics in 1989 from Princeton University 
and in 1995 the Ph.D. in physics from the University of California, San 
Diego. After a Department of Energy postdoctoral fellowship at the Oak 
Ridge National Laboratory, he joined the Institute for Nonlinear Science at 
UCSD as an Assistant Research Scientist. His research interests are in chaotic 
communication and the application of information theoretical algorithms to 
observed nonlinear dynamics. 



